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In this paper we review our recent work on the theoretical approach to quantum 
Loschmidt echoes, i.e. various properties of the so called echo dynamics - the composition of 
forward and backward time evolutions generated by two slightly different Ifamiltonians, such 
^ ■ as the state autocorrelation function (fidelity) and the purity of a reduced density matrix 

traced over a subsystem (purity fidelity). 



\ Our main theoretical result is a linear response formalism, expressing the fidelity and 

purity fidelity in terms of integrated time autocorrelation function of the generator of the 
perturbation. Surprisingly, this relation predicts that the decay of fidelity is the slower the 
faster the decay of correlations. In particular for a static (time-independent) perturbation, 
and for non-ergodic and non-mixing dynamics where asymptotic decay of correlations is 
absent, a qualitatively different and faster decay of fidelity is predicted on a time scale oc 1/5 
as opposed to mixing dynamics where the fidelity is found to decay exponentially on a time- 
scale (X 1/(5^, where 5 is a strength of perturbation. A detailed discussion of a semi-classical 
regime of small effective values of Planck constant h is given where classical correlation 
functions can be used to predict quantum fidelity decay. Note that the correct and intuitively 
expected classical stability behavior is recovered in the classical limit ft ^ 0, as the two limits 
\ 5^0 and — > do not commute. The theoretical results are demonstrated numerically 

^— «j . for two models, the quantized kicked top and the multi-level Jaynes Cummings model. Our 

■"".^ ' method can for example be applied to the stability analysis of quantum computation and 

' quantum information processing. 

P5 ■ §1. Introduction 

qh| In the early days of statistical mechanics, J. Loschmidt in a discussion with 

L. Boltzman suggested to study irreversibility by changing the velocities of all the 
molecules in a box which are, after an equal amount of time, supposed to return to 
^ ' their initial positions. Since the discovery of deterministic chaos, however, we know 

■ that such a gedanken experiment won't work after a short while since typical systems 

of many particles possess exponentical sensitivity on the variation of initial condition. 
Therefore arbitrarily small but non-vanishing perturbation of the trajectory (i.e. 
initial condition), or perturbation of the equations of motion (e.g. introducing a 
small external force field like gravity), drive the returning orbit away and hence the 
orbit will never return to the initial phase space point. Both mechanisms, namely 
perturbing the initial condition or perturbing the Hamiltonian, produce a similar 
effect in classical mechanics. 

This, however, is not the case in quantum mechanics. Due to manifest linearity 
and unitarity of quantum equation of motion (Schrodinger equation), quantum dy- 
namics is always stable against small variations of the initial state described by the 
wave-function.D Yet small variation in the Hamiltonian can produce interesting and 
highly non-trivial effects on quantum time evolution. 

It has been suggested by A. Fere^^ that susceptibility of quantum evolution 
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to small system variations, as measured by the Fidelity F{t) = of the 

states of the perturbed \ip5{t)) and unperturbed iV'(i)) time evolution, can provide 
a useful signature of classical chaos in quantum motion. Peres suggested that clas- 
sically chaotic systems are characterized by fast exponential decay of fidelity F{t), 
while the decay in regular systems should be qualitatively slower. This conclusion 
was asserted, based on a numerical experiment where an initial coherent state has 
been placed, respectively, inside the chaotic region of classical phase space, or in the 
middle of KAM island of classically regular motion. This view, elaborated semiclas- 
sically by Jalabert and Pastawski-- in the regime of initial coherent states and very 
short time scale, i.e. shorter than the so-called Ehrenfest time tE oc — log IS con- 
sistent with a semiclassical picture of decoherence of Zurekj^ which predicts, in the 
same semiclassical regime, that von Neumann entropy of a reduced density matrix of 
a central system traced over the environment grows with the rate proportional to the 
classical Lyapunov exponent (or better to say, local classical phase space stretching 
rate). According to this picture, fidelity decays exponentially F{t) tx exp(— At) for 
classically chaotic systems, with the perturbation independent rate A which matches 
the local classical phase-space stretching rate. 

We stress that this picture is justified under two rather severe assumptions: (i) 
coherent initial states which allow quantum-classical correspondence in phase space 
and (ii) short-times t < t-^ which guarantee pointwise quantum-classical correspon- 
dence of time evolution. If either of the two conditions is broken, the above picture 
can be shown to be incorrect.-' In quantum information processing, in particu- 
lar, one is certainly not interested in processing coherent initial states but rather 
random initial states which contain a maximal amount of quantum information. 
Furthermore, one is interested in the behaviour of fidelity F{t) for asymptotically 
long times, certainly longer than the Ehrenfest barrier ^e- 

In a series of paper^''^'^'^'^'^ we have elaborated on a linear response 
approach to fidelity decay which allows to deduce the behaviour of F(t) in the entire 
range of times for sufficiently small perturbation strength. In several cases linear 
response (perturbative) results can be extended to include all-orders in perturbation 
parameter and yield asymptotic long-time tails of F{t). The central result of this 
work is a linear-response fluctuation-dissipation formula which expresses fidelity F{t) 
in terms of integrated time-correlation function of the perturbation. As a general 
consequence of this formula, it follows that stronger correlation decay (typically as- 
sociated with stronger classical chaos of the underlying classical counterpart) means 
higher fidelity, or slower decay of fidelity. This has been quite an unexpected result 
as it seems just oposite to the short-time semiclassical picture.!^ However it can be 
shown, as a consequence of a delicate competition of time-scales (see e.g. Ref.^ or 
the present paper), that the two pictures are not contradictory but they're in fact 
complementary, and that there is a crossover of fidelity decay F{t) at t ~ from 
classical to quantum behaviour. 

Fidelity decay is sometimes associated with decoherence, which is a dynamical 
property of open quantum system coupled to another quantum system interpreted 
as environment. In order to elucidate this from a dynamical point of view, we 
have 

introducecPEDCl 

a novel concept of purity fidelity, namely the purity of a 
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reduced density matrix of a central system traced over the environment (or another 
part of a system) after undergoing the (Loschmidt) echo dynamics. We show an 
intimate relationship between purity fidelity and fidelity decay, and develop a linear 
response theory for purity fidelity in terms of special time-correlation functions of 
the generator of perturbation. 

This paper is a short comprehensive review of theoretical results on fidelity and 
purity fidelity which appeared in a series of recent papers.'^''^'^'^'^'^ Beyond 
that we discuss important aspects which have not yet been considered before, such 
as the case of coherent initial states, regular classical dynamics and 'ergodic' pertur- 
bation with vanishing time-average. 

It should be noted that a considerable recent interest in quantum Loschmidt 
echoeP'IHJ'inJmnniO,injniC10,IE),IT7JCH),El has been largely stimulated by spin 

echo experiments performed by the group of H. Pastawski,"^ and that fidelity is 
used as a benchmark for the quantum information processes.!^ 

§2. General theory of fidelity: linear response and beyond 

Let us consider a unitary operator U being either (i) a short-time propagator 
U = exp{—iHA/h) generated by some time-independent Hamiltonian H, or (ii) a 
Floquet map U = Texp(— i dTH{T)/h) of periodic time-dependent Hamiltonian 
H{t+p) = H{t). Arbitrary time evolution is generated by a group C/* where integer 
t is a discrete time variable. Note that in the continuous time case (i) we may let 
Z\ — > 0, so r = tA becomes a continuous time variable. Thus we shall consider the 
case of discrete time and develop general formalism for that case, whereas in the end 
the formulae for continuous time shall simply be obtained as the limit Z\ ^ 0. 

A general (static, for time dependent case see Ref.^) small perturbation of uni- 
tary propagator U can be 'parametrized' in terms of some bounded hermitian oper- 
ator V as 

U5 = Ueiq){-iV5/h), (2-1) 

where 5 is a small strength parameter. We note that V corresponds to the pertur- 
bation of the Hamiltonian = H + 5V in the continous time limit Z\ ^ 0. We also 
want to keep dependences on the (effective) Planck constant h explicit so we can 
control the semiclassical behaviour. For systems with a well defined semiclassical 
limit we shall be interested in the perturbations V which have well defined classical 
limit or Weyl symbol v{q,p). 

Starting from some initial reference state IV'), we consider two quantum time 
evolutions \ip{t)) = U^\ip), \'4'sit)) = U^\ip), and investigate a distance between the 
resulting Hilbert states as measured by fidelity F{t) = \{^ps{t)\tp{t))\'^ . It is crucial 
to note that fidelity can be written in terms of expectation value (•) = • {ip) of 
the unitary echo operator 

Ms{t) = Ui'U\ (2-2) 



namely 



F{t) = \{Msm'- 



(2-3) 
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We shall later discuss the dependence of fidelity on dynamical properties of the time- 
evolution U, on the strength of perturbation 5, and on the structure of initial state 

Concerning the last point it is often useful to study fidelity with respect to 
a random initial state \^p), i.e. to average, denoted by ((•)), over an ensemble of 
initial states with a unitarily invariant measure in A^-dimensional Hilbert space. 
This should generally correspond to processing of states with maximal quantum 
information (entropy) which is certainly the situation of most interest for quantum 
computation.^^' We note that bilinear expressions in 1-0) are averaged by means of 
a trace, 

mA\i^))^=--m = ^tTA, (2-4) 

however higher order expressions, like fidelity (|2-3() which is of fourth order in random 
variable have to be averaged by means of pair contractions in the asymptotic 
Gaussian (M — > oo) limit. For example, one can compute the average fidelity {{F{t))) 
in terms of the average fidelity-amplitude {{f{t))) = {{Ms{t))) with a small, semiclas- 
sically vanishing correction 

{{Fm = \{{fmf + ^. (2-5) 

Thus we have shown that the variation of fidelity with respect to random initial states 
becomes unimportant for large Hilbert space dimensions M, e.g. when approaching 
either the semiclassical or the thermodynamic limit, so one is justified to average the 
amplitude f{t) in order to simplify theory.'^ 

On the other hand, for the purpose of quantum classical correspondence, it is 
most suitable to consider coherent initial states {ip), in some sense the states of 
^minimal quantum information\ Below we shall specify our general formulae for the 
two extremal cases, namely random and coherent initial states. 

To derive our general theory we first rewrite the echo operator 1)2 -2(1 in terms of 
a time-dependent perturbation operator in the interaction picture 

Vt := U-^VUK (2-6) 

Then we recursively insert the expression of unity 1 = f7* in the definition 
(tr^ and observe U~^'U^^U^+^' = exp{iVtf5/h) for t' running from t — 1 downto 0, 
yielding 

= C/~^*"^^?7*-iexp(il4_i(5/;i) = 

= U^'^'''^'^U'~^e^v{iVt-25/h)Qx:p{iVt-i5/h) = 

= e^^{iVo5/h) e-x.^{iVi5/h) ■ ■ ■ eiq){iVt-i5 / K) . (2-7) 
The obvious next step is to expand the product (|2-7|) into a power-series in 5 

in=l tl,...,tm=0 
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where the operator T denotes a left-to-right time ordering. Such a perturbative 
expansion converges absolutely for any S provided that the perturbation y is a 
bounded operator. Therefore fidelity may be computed to arbitrary order in 5 by 
truncation of the above expansion (|2-8() and plugging it into ()2-3() . So we can see 
that the fidelity F{t) can be expressed entirely in terms of multiple time correlation 
functions of the generator V of the perturbation. 

To second order in S one obtains a very useful linear response formula 

^(*) = l-5 E C{t',t") + 0{6'), (2-9) 

t',t"=0 

where 

C{t',t") := {Vt^Vt,,) - {Vt,){Vt,,) (2-10) 

is a 2-point time correlation function of the quantum observable V. This formula can 
be interpreted in terms of a dissipation-fluctuation relationship. On the LHS we have 
fidelity which describes dissipation of quantum information and on the RHS we have 
an integrated time-correlation function (fluctuation). A simple-minded qualitative 
conclusion drawn from the formula (|2-9|1 says: The stronger the decay of correlations 
the slower the decay of fidelity and vice versa. The linear response formula can be 
rewritten in a slightly more compact notation as 

F{t) = 1 - ^CintW, (2-11) 

where 

t'=0 

We shall use this compact notation later in section |31 

Note that the range of validity of the above formula (|2-9I2-11|) is not at all re- 
stricted to short times. The only condition is that 6 is sufficiently small such that 
1 — F{t) is small (<C 1). Below we shall discuss different regimes and different time- 
scales based on this formula and its higher order extensions. In the semiclassical 
regime of approaching the classical limit /i — > the quantum correlation function 
C{t',t") goes over to the corresponding correlation function of the classical pertur- 
bation v{q{t),p{t)). 

Therefore the fidelity for classically chaotic systems will decay with the rate 
which is inversely proportional to their rate of mixing. Furthermore for classically 
non-ergodic, i.e. regular or integrable motion, the correlation functions will generally 
not decay to zero and the fidelity will therefore decay much faster. 

For continuous time t one expects that in the so-called Zeno regime of short 
times, such that correlation function does not yet decay appreciably. Fidelity will 
always decay quadratically F{t) ~ 1 — (5^C(0)t^//i^. This is independent of the nature 
of the corresponding classical dynamics, whether it is regular or chaotic. However 
we are interested in longer times, beyond the range of quantum Zeno dynamics. 
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2.1. Regime of ergodicity and fast mixing 

Here we assume that the system is (classically) ergodic and mixing such that 
the correlation function C{t',t") decays sufficiently fast as \t' — t"\ grows; this typ- 
ically corresponds to globally chaotic classical motion. This also implies that after 
a certain short (Ehrenfest) time scale t^, which can for a chaotic system be written 
~ ^og{l/h)/X in terms of an effective classical Lyapunov exponent A, the time 
correlation functions become independent of the initial state and thus equal to the 
random state average, e.g. 

Cit, t + T) = C{t) := {{VVr)), if t > tE. (2-13) 

In this situation we can safely assume random initial states, or argue that for times 
larger than tE the results are the same for any inital state. Then the linear response 
formula can be rewritten as 

= 1-^1 2*^(0) + E - *')^(*') I + ^('^')- (2-14) 

Now we shall assume that correlation function C{t) decays sufficiently fast, i.e. faster 
than 0(t~^) and that a certain effective time-scale tmix of decay of C{t) exist. For 
times t ^ tmix we can neglect the second term under the summation in (|2.14() and 
obtain a linear decay in time t in the linear response regime 

1 °° 

Fem{t) = 1 - 2{5/hfat with a = -C(0) + E C{t). (2-15) 

t=i 

Here, a plays the role of a transport coefficient. 

We can make a stronger statement valid beyond the linear response regime if we 
make an additional assumption on the factorization of higher order time-correlations, 
namely that of n— point mixing. This implies that 2m-point correlation {{Vt^ ■ ■ ■ Vtj^)) 
is appreciably different from zero for t2m — ti — > oo only if all the (ordered) time 
indices {tj,j = 1 . . . 2m} are paired with the time differences within each pair t2j — 
t2j-i being of the order or less than tmix. Here we have assumed without loss of 
generality that the perturbation is traceless {{V)) = 0. If not, we subtract {{V))l from 
V which does not affect probability F{t). Then we can make a further reduction, 
namely if t 3> mtmix 

T E {{yt.Vt,---VU)^f E {{VM)---{{Vt,^_,Vt,J)^^^{tar. (2-16) 

tl,...,t2m=0 tl,...,t2m=0 

We then obtain a global exponential decay 

i^cm(i) = \{{Ms{t)))f = exp {-2t/T,^), Ten, = (2-17) 

with a time-scale rem = 0{6~'^). We should stress that the above result (|2-17j) has 
been derived under the assumption of true quantum mixing,'23 which can be justified 
only in the limit AA — > oo, e.g. either in semiclassical or thermodynamic limit where 
the correction H2-5|) is irrelevant. Note that the result H2-17() can be connected to 
Fermi golden rule^ as it is based on time-dependent perturbation theory. 
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2.2. Non-mixing and non-ergodic regime 

The opposite situation of non-mixing and non-ergodic quantum dynamics, which 
typicahy corresponds to integrable, near-integrable (KAM), or mixed classical dy- 
namics, is characterized by a non-vanishing time-average of the correlation function 

^=;i™^ E ^(*''*")- (2-18) 

t',t"=o 

Here, due to non-ergodicity, the time-average C depends on the structure of the 
initial state {ip). We further assume that a certain characteristic averaging time- 
scale tave exists, namely it is an effective time t = tave at which the limiting process 
H2-18() converges. Therefore, for sufficiently large times t » iavej the double sum on 
RHS of eq. (|2-9() can be approximated as Ct'^, so the linear-response formula H2-9() 
yields, in contrast to H2-15() . a quadratic decay in time 

F.e{t) = l-( — y + 0{6% ^ne = -^, (2-19) 

with time-scale r^e = 0{6^^). One should observe that the non-ergodic time-scale 
Tne can be much smaller than the ergodic -mixing tiiii6-sca.l6 Tq^yi 

(|2-17|) provided h is 

fixed, or the limit (5 — > is taken prior to the limit ^ ^ 0. Yet Tne is typically still 
much longer than the Zeno time scale of universal quadratic decay. 

Again we can make a much stronger general statement going beyond the second 
order 5-expansion. If we assume that t ^> tave, we can re- write the m-tuple sums in 
the series ()2-8|) in terms of a time average perturbation operator 

t-i 

V = lim (l/t)y 14/, (2-20) 

t'=0 

namely 

= E 1^^^") = ^^^P (^^^V^)>, F^eit) = \fneit)\'. (2-21) 
m=0 

Note that V is by construction an integral of motion,'^'^^ \U-,V\ = 0, and reduces 
to a trivial multiple of identity in the case of ergodic dynamics studied in previous 
subsection. Whereas in an ergodic and mixing case, m— th order term of (|2-8p grows 
with time only as 0{t"^^'^) (for even m) since it is dominated by pair time correlations, 
here in a non-ergodic case, the non-trivial time average operator V already gives 
the dominant effect, namely 0{t'^) for m— th order term of (|2-8|) . so the effect of 
pair time correlations can safely be neglected for sufficiently long times (i 3> iave)- 
Observe also that time averaged correlation is just a variation of the time averaged 
perturbation 

C=(F2)_(y)2_ (2-22) 
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In the semiclassical regime for random initial states, C goes to a purely classical 
(^^-independent) quantity Cd = {v'^)ci — where is a time-averaged classical 
limit of observable V and (•)ci is a classical (microcanonical) phase-space average. 

Our conclusions may not be valid in the special case where the time-average V 
is a trivial operator, namely when it either vanishes or is proportional to identity. 
This can happen for very special choices of perturbations V or for systems and 
perturbations with particular geometric or algebraic symmetries. This option also 
implies that the classical time-average is trivial v = const, and that (7 = 0. Of course, 
in such a case, fidelity decay has to be discussed separately, see e.g. Refs.^^''^ 

Let us now use expression H2-21() to derive some explicit semiclassical results 
in the special case of integrable classical dynamics. For a system with d degrees 
of freedom we thus have d canonical constants of motion - the action variables 
/ = (/i, . . . , Id), which are quantized using EBK rule In = h{n -|- 7/4) where n = 
(ni, . . . , rid) is a vector of integer quanum numbers rij and 7 = (71, ... , 7^) is a vector 
of integer Maslov indices 7^-; the latter are irrelevant for the discussion that follows. 
Since the time averaged operator V commutes with U and with the actions /, it is 
diagonal in the (generically non-degenerate) basis of eigenstates of / (quantized tori) 
|n). In leading semiclassical order one may write 



where v{I) is the corresponding classical time-averaged observable in action space. 
The fidelity H2-21|) can therefore be written as 



Provided the diagonal elements of the density matrix can be written in terms of 
some smooth structure function D{In) = \ {n\ip)\'^, and replacing the sum (|2-24|) by 
an integral over the action space, which is justified for small h up to classically long 
time oc hP6~^, we obtain 



The obvious next step is to compute this integral by a method of stationary phase. 
However the result depends on the precise form of the function D(I) which may in 
turn depend explicitly on h. Below we work out the details for two important special 
cases, namely a random and a coherent initial state. 

2.2.1. Semiclassical asymptotics for a random initial state. 

Let us first assume uniform averaging over (random) initial states D(I) = 1/M = 
{2-Kh)'^ /V. For large t5/h the above integral H2-25() can be written as a sum of contri- 



butions stemming from, say p points, I-q^r] = 1, . . . ,p where the phase is stationary, 
dv{lr,)/dl = 0. This yields 




(2-23) 




(2-24) 



n 




(2-25) 



fZit) 



(27r)3'^/^ n_ "^/^ A exp{itv{I^)6/h + iiyr^} 
V t6 ^ |detVJV2 



r,=l 



(2-26) 
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where {V^jj-jt := d'^v{Ir^) / dijdik is a matrix of second derivatives at the stationary 
point r/, and = 'ir{m-^ — m_)/4 where m± are the numbers of positive/negative 
eigenvalues of the matrix V^. The stationary phase formula (|2-26|) is expected to be 
correct in the range const h/5 < t < const'/ J. Most interesting to note is the asymp- 
totic power-law time and perturbation dependence -Fn™ ~ which allows 
for a possible crossover to a Gaussian decajJ^ when approaching the thermodynamic 
limit (i — > oo. 

2.2.2. Semiclassical asymptotics for a coherent initial state. 

Now let us consider a single d-dimensional general coherent state centered at 
{I*, 6*) in action-angle space 

{n\r,e*) = (^^y\detA\'/' exp - I*) • ^(/^ _ /*) - m • r | , (2-27) 

where yl is a positive symmetric d x d matrix of squeezing parameters, giving 

D{I) = {h/7ry/^\detA\^/^eM-{I-I*)-Ml-I*)/f^) (2-28) 

and 

|detyl|^/2 



jd'^I exp !^^{itv{I)6 - {I - I*) ■ A{I - 1*))^ . (2-29) 



Using the assumption 6t <C 1, we see that a unique stationary point Is of the 
exponent approaches /* as 5 — > 0, 

I^ = r -!^A-^v' + 0{5''), where v' := (2-30) 



so we may explicitly evaluate H2-29() by the method of stationary phase without any 
lower bound on the range of time t, 

Note that the fidelity decay for a coherent initial state with regular classical motion 
has a time-scale 



''"ne— coh 



{2h)'/\v' ■ A-W)-^/^5-' oc h^/^6-\ (2-32) 



which is consistent with (|2-19() with C = ^h{v' -A^^v') and is by a factor proportional 
to h^^/'^ longer than the time-scale t^q (x h/5 for a random initial state. It should 
be noted that the above derivation of fidelity decay for a coherent initial state l\'2-27\ - 



12 -321) remains valid in a near-integrable (KAM) situation of mixed classical phase 
space, provided that the initial wave packet is launched in a regular region of phase 
space where (local) action-angle variables exist. 
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2.3. Finite size effects and time and perturbation scales 

The theoretical relations of the previous subsections are strictly justified in the 
asymptotic limit M ^ oo. For finite M, fidelity F(t) cannot decay indefinitely but 
starts to fluctuate for long times due to discreteness of the spectrum of the evolution 
operator U. Let us write the eigenphases of U and Us, and the corresponding 
eigenvectors, respectively, as 4>n^ and \(l)n), \4'n)^ n = 1, . . . ,Af, satisfying U\4>n) = 
e~*'^"|(/)„), Us\(j)n) = e~*'^"|(/)^). Now define a unitary operator W which maps the 
eigenbasis of U to the eigenbasis of Us, namely := \(f)f^) for all n, with matrix 

elements Wmn ■= {4'm\W\(j)n) , and write the initial state as ipn = {(pn\tp)- Note 
that the matrix Wmn is real orthogonal if U and Us possess a common anti-unitary 
symmetry (e.g. time-reversal). Now it is straightforward to rewrite the fidelity 
amplitude as 

fit) = E ^ni^MWmi exp (i(</>„ - cAf )t) , (2-33) 

l,m,n 

At this point we are interested in the long time fluctuations so we compute the time 
averaged fidelity fluctuation 

Ft. := \hn ^ r dt\f{t)\^= T rkW^i^mWmiMWnil''. (2-34) 

k,l,m,n 

In the process of averaging over the time we have assumed that the eigenphases 
are non-degenerate so exp{i{(f>n — (pn' + 4*1 ~ 4>f')t) = ^nn'^i,!'- We see that fidelity 
fluctuation Fta depends on the orthogonal/unitary matrix Wmn and the initial state. 
Detailed discussion of various cases are given in Here we only give results 

for random intitial state, where ipn can be assumed, for large M, to be independent 
complex random Gaussian variables with variance 1 /M. Averaging over an ensemble 
of initial states yields 

{{Ft.)) = ^\ Y1 1^"'!' + Yl \Wni\''\Wmi\'' I . (2-35) 



nJ n,ml 



For sufficiently weak perturbation, one may assume that the matrix Wn,m is near 
identity Wn,m ~ ^n,m, thus -F'ta-^'f ~ 2/AA while for strong perturbation Wn,m can 
be assumed to be a random matrix, yielding Fl^°^^ = l/M. Therefore the l/M 
dependence of the residual fidelity Ft., for random intial states does not depend on 
the structure of eigenstates W^^m^ and there is only a factor of 2 difference between 
the two extreme cases. 

In addition to a random initial state, we shall now assume ergodic and mixing 
classical dynamics, so fidelity decay is initially given by exponential law (|2-17|) until 
it reaches the plateau (|2-35() . The saturation time t^ where -Fem(i*) = -^ta, is the first 
important time scale 

1 fPd 
t, = -T^^lnM^j^H^/h), (2-36) 

where ad is the classical limit of the transport coefficient (|2-15|) . 
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The second new time-scale is related to the asymptotic non-decay of time correla- 
tions for finite- quantum dynamics, namely even if the system is classically mixing 
the quantum correlation function will have a small non- vanishing (?i-dependent) time 
average 

C = WWW)) = E \^n\\Vnn)^ (2-37) 
n 

where Vmn = {fprnlylfpn) ■ Since, however, the classical system is assumed to be 
chaotic, implying ergodicity and mixing, the matrix elements V^„ behave like pseudo- 
random variables with a variance given by the Fourier transformation S{uj) of the 
corresponding classical correlation function Cci(i) at frequency uj = (j)m ~ 4'n- 

On 

the diagonal we have u; = and an additional factor of 2 due to the form of the 
invariant random matrix measure (see e.gP^). Thus we have 

C = ^, (2-38) 

where iTcI = 5'(0)/2 is the classical limit of (|2-15() . The result, due to ergodicity, 
does not depend on the particular form of the initial state The decay of fidelity 
H2-14() will start to be dominated by the average plateau H2-38() at sufficiently long 
time t, when X^^/I^^g ~ ^')^ ~ 2<Jc\t^ /M > (Tdt, i.e. for times t greater than tp 

= IaA oc h-'^, (2-39) 

which is just the Heisenberg time associated to the inverse density of states. 

Depending on the interrelation among four (or five) time-scales r^m oc h?5~'^^ 
tp oc h'''^6^, oc h^ln{l/h)6^'^d, tmix h^S^, (and t-^ oc ln(l//i)(^° if we are con- 
sidering coherent initial states, like e.g.'^'^'^) we can have four (or five) different 
regimes depending on the three main scaling parameters: perturbation strength 6, 
Planck's constant h, and dimensionality d. Note that we always have t^, > rem- All 
regimes can be reached by changing only the parameter 6 while keeping h and d fixed 
(see fig. CJ: 

(a) For sufficiently small perturbation 6 we will have tp < rem- This means that 
Fem(tp) is still close to 1 and we will initially have quadratic decay (|2-19jl with C 
given by an average finite size plateau 1)2-38(1 . This will occur for 5 < 6p where 

6 -n(^\'" -^^^hfl--^ (2-40) 

In fact, in this regime, also referred tc^^'EU as perturbative, one may use first order 
stationary perturbation theory on the eigenstates of Us, yielding (/)^ = + 14„5//i+ 
0(^2), Vmn = Smn + 0{6), and rewrite (following) the finite size fidelity (|2-33() in 
terms of a Fourier transform of a probability distribution w{V) of diagonal matrix 
elements Vnn, /pert(i) = / dVw{V) eyip{—iVt6/h). Since wiV) is conjectured to be 
Gaussian for classically ergodic and mixing system,!^ it follows that /pert(i) is also 
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a Gaussian with a semiclassically long time-scale Tp 

/pert(i)=exp(-(t/Tp)V2), rp=^- =^_^^^. (241) 



(27r)^/2c 

(b) If < fp < t* we will have a crossover from initial exponential decay of fidelity 
(|2-17j) to a Gaussian decay (|2-4H) at t ~ tp, which will terminate and go over to 
fluctuating behavior when -Fpcrt(i) = -Pta(AA). Note that this will happen before 



time t=K which is estimated based on a slower exponential decay (|2-17|) . This regime 
will exist in perturbation range 6p < 6 < 6s with an upper border 6s determined by 
the condition tp = t* to be 

<5s = (lnAA/2)i/2^p. (2-42) 

(c) If we still further increase 6, we have the most interesting, 'fully nonperturbative' 
regime, when < t=K < tp and we will have a full exponential decay H2-17() . up to 
time t=K when the fidelity reaches finite size fluctuations. This regime continues as 
long as 6 < (5mix where the border 

* ' ^ (2.43) 



V Ccl^mix \2t 



•rnix 



is determined by the condition rem = ^mix which is a point where the arguments 
leading to the factorization H2-16() and exponential decay (|2-17|) are no longer valid. 
We note that the relative size of this window range (5mix/5s = -y/AA/ (tmix InAA) in- 
creases, both, in the semiclassical and in the thermodynamic limit. This regime also 
corresponds to 'Fermi golden rule decay' discussed in.^^ 

(d) Further increasing 6 > ^mix, the estimated fidelity decay time eventually be- 
comes smaller than the classical mixing time t^ix; '^em ^ ^mix* this regime, the 
perturbation is simply so strong that the fidelity effectively decays within the shortest 
observable time-scale (tmix). 

However if we consider a non-random, e.g. coherent initial state then the quan- 
tum correlation function relaxes on a slightly longer, namely Eherenfest time-scale 
tE so the regime (c) should terminate already at a little smaller upper border 6 = 6e 
which is naturally determined by the condition Tem = 

^[aeiln(l/n)]i/2 ~ [ln(l7r)]V2- (2-44) 

For coherent initial states one thus obtains an extra but very narrow regime 6-e < 
S < (5mix (describing the time-range tmix < t < ts) where the fidelity decay can be 
computed in terms of classical Lyapunov exponents.^''^^' 

Yet in the regime of non-ergodic, say integrable classical mechanics things are 
simpler, as we do not have to worry about the average plateau in the correlation 
function due to a finite M because we already have a higher average time correlation 
C Cci 7^ (|2-18|) . Thus we have here only two relevant time-scales, namely Tne 
giving initial quadratic decay ()2-19|) . and the saturation time-scale t^,, due to finite 
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Fig. 1. Schematic view of different regimes (a-d) of fidelity decay in the regime of classically chaotic 
(or mixing) dynamics. 



size fidelity fluctuation, which depends on the properties of the initial state (power 
law ()2-26() for a random initial state, versus Gaussian H2-31() for a coherent initial 
state). However in the particular special case of 'ergodic' perturbation V = oi 
regular dynamics (as mentioned above), one now clearly obtains that (7 = which 
is equivalent of having all diagonal matrix elements vanishing, i.e. Vnn = [see eq. 
H2-37() ]. In this case the time scale of fidelity decay Tne from eq. H2-19() formally 
diverges. This in fact means that fidelity decays much slower than in generic case, 
now the decay being described by the first nonvanishing - i. e. fourth order term in 
the ^-expansion.l^ 

Note that our result is not contradicting any of the known facts of the quantum- 
classical correspondence. For example, a growth of quantum dynamical entropie^^ 
persists only up to logarithmically short Ehrenfest time which is the upper bound 
for the validity of the strict classical-quantum correspondence in the fidelity decay as 
derived in Ref.^''^'*' and within which one would always find F^°^{t) > Fem{t) above 
the perturbative border 5 > 6p, whereas our theory reveals new nontrivial quantum 
phenomena with a semiclassical prediction (but not correspondence!) much beyond 
that time. On one hand, if we let ^ ^ first, and then 5 — > 0, we recover a result 
supported by a classical intuition, namely that the regular (non-ergodic) dynamics 
is more stable than the chaotic (ergodic and mixing) dynamics. On the other hand, 
if we let (5 ^ first, and only after that ?i ^ 0, we find somewhat counterintuitive 
results saying that chaotic (mixing) dynamics is more stable than the regular one. 



§3. Purity fidelity 



In this section we propose another characteristic which describes the quality of 
Loschmidt echoes, in particular for systems composed of two parts corresponding to 
different degrees of freedom, for example the central system and the environment. 
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We consider a system with a product Hilbert space 7i = TCi "^2, where 7ii^2 
are factor Hilbert spaces of the two subsystems. We are again studying an echo 
experiment, now with the the initial state beeing a product state 

IV'(O)) = l^i) iV'a), (3-1) 
yielding a, generally entangled, final state 

\m) = MsitMm- (3-2) 

The entanglement produced by imperfect quantum echo can be most directly char- 
acterized by the reduced density matrix 

pi{t) = tT2\m){m\ = tr2 (M5(t)|^(0))(V(0)|M](t)) (3-3) 

obtained by tracing out the second subsytem 7^2- Now one can say that iV'(i)) is 
disentangled iff the density matrix pi (t) is pure, and entanglement can be character- 
ized by purity; specifically for the echo situation we call this quantity purity fidelity 

Fp{t)= tvi[pi{t)f. (3-4) 

Note that purity fidelity is in some sense a weaker quantity than fidelity: in order 
for fidelity to be high, purity fidelity must also be high, but not vice versa. High 
purity fidelity only requires the state to return to the factorized form, which is a 
much weaker condition than complete recurrence as required by high fidelity. 

Now we apply our perturbative expansion of the echo operator H2-8() to the 
expression for purity fidelity (|3-4() . Keeping only terms up to second order in the 
perturbation strength 6 we obtain a linear response formula for purity fidelity, 

Fp{t) = 1 - 26^h-^ {Ci,t(t) - Ant(t)} + • • • , (3-5) 

Ant(t) := 5] 1(1, i^imih i)p + i\m\i, 

which is analogous to a simpler formula (|2-llj) for fidelity. We have used an explicit 
notation of a complete basis of a product Hilbert space as \i, v) = \i) where latin 
indeces label states in the first and greek indeces states in the second subspace. There 
is a general relation between fidelity and purity fidelity which can be formulated in 
terms of a rigorous mathematical inequalitjJ^ 

Fp{t) > F\t). (3-6) 

The formal similarity of linear response formulae (|2-llj) and (|3-5|) results in 
similar physical behaviour of fidelity and purity fidelity, in both qualitatively different 
cases of dynamics, ergodic and mixing — chaotic, and regular. For example, in the 
regime of ergodic and mixing dynamics we can use the result that fidelity is after a 
while, due to ergodicity, independent of the initial state, hence we can write for the 
echo operator in the weak limit sense 



Ms{t) ^ exp{-t/Tem)l. 



(3-7) 
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Plugging this expression into formula (|3-3|) and the obtained result into formula 
H3-4(l we finally obtain a simple expression compatible with the upper bound of the 
inequality 1)3 -Ol) 

Fp(f) =exp(-4t/rem). (3-8) 
In this case the relative weight of the correction term Dint in the linear response 
formula (|3-5|1 vanishes in proportion to Cjnt in the semiclassical limit /i — > 0, — > oo. 
Again, the same consideration of time and perturbation scales applies for finite M 
as discussed above for the case of fidelity. 

In the non-mixing or classically regular case, the situation is more complicated 
as purity fidelity depends on the structure of initial state The only thing we can 
state generally in this case is the quadratic decay in the linear response regime, i.e. 
as long as 1 — Fp(t) is small we have 

Fp{t) = l--^{C-D)t^ + 0{t^). (3-9) 

Since in this case time-averaged correlation functions are non-vanishing, the inte- 
grated ones grow as oc t^, so we have written Cint(i) = Ct'^, Diat{t) = Dt^. It is 
interesting to discuss the two extreme cases of intial states: (i) In the case of coher- 
ent initial states (Gaussian wave packets) we have showrPJ that the term D cancels 
the term C in the leading semiclassical order, so C — D cc h? , whereas C (x h. This 
means that purity fidelity for coherent initial states decays with an h independent 
time scale which is by a factor proportional to h~^/'^ longer than the time scale of 
fidelity decay, (ii) On the other hand, in the case of initial random states one can 
show that the term D is negligible compared to C in the semiclassical limit, so both, 
fidelity, and purtity fidelity decay with the same rate, similarly as in the case of 
chaotic dynamics. 

There is an important special case of potential practical interest, if the perturba- 
tion is such that the unperturbed composite system is decoupled^ i.e. Us=q = Ui®U2. 
Then the unperturbed (e.g. forward) evolution does not change the purity, so the pu- 
rity fidelity (purity of echo dynamics) is the same as purity (or linear entropy) of uni- 
directed time-evolution alone. Again, our linear response formalism predicts faster 
increase of linear entropy (decay of purity) for regular or weakly chaotic systems 
than for strongly chaotic systems.^ This result has been independently reproduced 
in Ref.123 

In next section we outline some of our numerical results which confirm the theory 
of the last two sections. 

§4. Numerical experiments 

We shall consider two numerical toy models by which we may demonstrate and 
verify the theoretical results of previous sections. 

First we choose Haake's quantized kicked topP^ since this model served as a 
model example for many related studies, see e.g. Refs.E3'^''^''22.i'J3i.34.i;35i 

unitary propagator of the kicked top reads 

U = U{a,-f) = exp{-ijJy)exp{-iaJ^/2J), (4-1) 
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where Jk {k = x, y, z) are quantum angular momentum operators obeying [J^, Ji] = 
i^kirJr- The (half)integer J determines the size of the Hilbert space 2 J + 1 and 
the value of the effective Planck constant h = 1 /J. The perturbation is defined by 
varying the parameter a, Us = U{a + S,j), so V reads 



The classical limit is obtained by letting J = 1/h ^ oo and writing the classical 
angular momentum in terms of a unit vector on a sphere r = {x,y,z) = J/ J. The 
Heisenberg equation for the SU(2) operators J, J' = WJU, reduces to the classical 
area preserving map of a sphere 



For a = the system is integrable, while with increasing a there is a transition 
to chaotic motion. The second parameter 7 is usually set to 7r/2, however in our 
numerical simulation we will use two different values exhibiting qualitatively dif- 
ferent correlation decay (for large a): the 'standard' case a = 30,7 = where 
Cci{t) decays in oscillatory way and the case a = 30,7 = where Cci{t) decays 
monotonically (see fig. 0)) . 

In the case of 7 = 7r/2 we have two discrete symmetries. The evolution U 
commutes with and Ry, the rotations of vr around the x and y axes, respectively. 
The Hilbert space is therefore reducible into three invariant subspaces (using notation 
of Peres's boolP with the basis \m) of eigenstates of Jz and assuming J to be an even 
integer): EE of dimension J/2 + 1 with the basis states |0) and {|2m) + 1 — 2m)}/\/2; 
00 of dimension J/2 with the basis {|2m — 1) — | — (2m — l))}/\/2; OE of dimension 
J with the basis {|2m) - | - 2m)}/^/2 and {|2m - 1) + | - (2m - l))}/\/2 with 
m = 1, . . . , J/2 in all three cases. For 7 7^ 7r/2 the spaces 00 and EE coalesce as 
Ry is the only discrete symmetry left. In numerical experiments we always choose 
the OE subspace so that the dimension of the Hilbert space is M = J. 

We will compute the fidelity of two different types of initial states: (1) ran- 
dom initial state with components Cm = ("ilV') being independent Gaussian pseudo- 
random numbers, and (2) pure minimal wavepacket initial state, namely SU(2) co- 
herent state = \'&,^) centered at the point n = (sin -i? cos sin -i? sin cos ??) on 
a unit sphere 




(4-2) 



x' = cos 7(x cos az — y sin az) + z sin 7 

y = y cos az + X sin az 

z' = z COS 7 + sin 7(7/ sin az — x cos az) . 



(4.3) 



Note that in the classical limit the perturbation generator is 



(r) = — . 

^ ^ 2 



(4-4) 




1/2 



cos {'d/2) 



J+m 



sin (??/2) 



J—m 



(4-5) 
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As a second model we choose a Jaynes-Cummings (JC) Hamiltonian, a pop- 
ular model in the realm of quantum optics. JC model is an autonoumous (time- 
independent) system of an harmonic oscillator interacting with a rotor, e.g. one 
mode of electromagnetic field, and a spin, say of an atom. Suppose the oscilla- 
tor is described by standard anihilation/creation operators a,a^, and the spin of 
(half)integer value J by SU(2) variables J±, Jz- Then the JC Hamiltonian reads 

H = hioa^a + heJz + (G{aJ+ + J_) + G'(aJ_ + J+fj . (4-6) 

Here the last (counter rotating) term has been included to allow for chaotic mo- 
tioiPl Contrary to the kicked top model and the framework of previous theoretical 
discussion, time t is here a continuous variable. There is, however, a natural limiting 
procedure in order to obtain continuous time versions of results of sections 2 and 3 
by replacing all the sums like X^^/I^o integrals dt' . A convenient scaling towards 
the classical limit J — > oo, ^ 0, is obtained by fixing classical angular momentum 
hJ = 1, hence h = 1/J. We perturb JC model by detuning., i.e. slightly changing 
the magnetic field parameter e, so the generator of perturbation reads V = hJ^. As 
an initial state we here choose only the coherent state, namely a direct product of 
oscillator coherent state with complex parameter a, and SU(2) coherent state with 
another complex parameter r = e*"^ tan(T?/2): 

IV') = e"'^'-"*"|0)2®(l + |T|2)-V^iJ,J)i. (4-7) 

The classical JC Hamiltonian is integrable if either G = or C = 0, whereas it can 
display all the variety of KAM regimes as both parameters G, G' are increased and 
reaching practical ergodicity for sufficiently large G,G'. We study numerically two 
different sets of parameters: In both we choose lo = e = 0.3, and (i) G = 1,G' = 
corresponding to classically regular motion, and (ii) G = G' = 1 corresponding 
to classically almost fully chaotic motion. In both cases we choose the same initial 
coherent state with parameters a = 1.15, r = 0.29-1- iO. 46 and set angular momentum 
J = 4. We first choose a small perturbation parameter 6 = 0.005. Numerical results 
of time correlation functions and fidelity decay in such a linear response regime are 
shown in figs. I2I3I 

In order to reach beyond the linear response regime and go deeper into the 
semiclassical regime we use the kicked top model. There we may hope to exploit 
the behaviour of classical correlation functions in order to predict quantum fidelity 
decay. First, we turn to the regime of completely chaotic classical dynamics. 

The classical correlation functions calculated by using the classical map (|4-3() 
are shown in fig. |1J For a = 30,7 = '^/'^ the correlation function is oscillating 
with an exponential envelope hence the transport coefficient ad = 0.00385 is quite 
small. For a = 30,7 = t^/^ the correlation decay is monotonic and exponential with 
(Tci = 0.0515. The decay of quantum fidelity 1)2-17^ can now be obtained by using 
the classical limit a ^ cr^. 



Fem(t) = exp(- 2,5 Welt). 



(4-8) 
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Fig. 2. Integrated quantum correlation function of the perturbation in the JC model, divided 
by time t. Solid curve corresponds to classically chaotic case (see text) and is approaching 
a constant plateau (dashed), whereas dotted curve corresponds to integrable case (see text) 
following linear growth (dashed). 




Fig. 3. Fidelity decay of the coherent initial state in JC model for regular and almost fully chaotic 
dynamics (see text). Dotted and full curves correspond to numerically exact evaluations of 
fidelity, whereas circles and squares reproduce a precise evaluation of the linear response formula 
H2-9|l in terms of quantum correlation functions (of fig.EJ. Dashed curves are best fitting linear 
and quadratic decays. 
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Fig. 4. Classical correlation function Cci(t) for kicked top with a = 30, and 'y — n/6 (top, full 
curve) and 7 = 7r/2 (bottom, broken curve). Finite time integrated correlation function is 
shown in the inset, converging to a^i = 0.00385 and 0.0515, for 7 — 7r/2 and n/6, respectively. 
Averaging over 10^ random initial conditions on a sphere is performed. 



This formula has been compared with the exact numerical calculation of fidelity 
averaging over a set of random initial states. We stress that we found no observable 
difference for sufficiently large J when we have instead chosen a fixed coherent initial 
state, however this calculation is not shown in the figures. As the finite size fidelity 
fluctuation level -Fta decreases with increasing Hilbert space dimension, we chose 
large J = 4000 in order to be able to check exponential decay (|4-8() over as many 
orders of magnitude as possible. The results are shown in fig. [21 The smallest and 
the largest 5 shown, roughly correspond to borders 5^ and ^mix, respectively. As 
we can see, the agreement with an exponential decay is excellent, at least over four 
decades in the fidelity F{t). Note that for 7 = 7r/2 and the largest 5 = 2> ■ 10~^ the 
time-scale of the decay of fidelity is comparable to the classical time-scale tmix so 
the factorization assumption (|2-16p is strictly no longer applicable. However, due to 
oscillatory correlation decay, overall agreement with the theory l|4-8|) is still rather 
good, but the oscillations of the correlation decay are reflected by oscillations of the 
fidelity decay (around the theoretical exponential curve). Of course one does not 
need such a large J in order to have an exponential decay, but for smaller J the 
fluctuation level Fta will be higher so the exponential decay H4-8|) will persist for 
correspondingly smaller time. 

Then we focus on the so-called perturbative regime 5 < 5p where the fidelity 
decay will be dictated by a finite size correlation average H2-38() . so according to eq. 
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Fig. 5. Quantum fidelity for a kicked top with parameters a — 30, J = 4000 and full trace average. 
Top figure is for 7 = tv/2 and for S = 5 ■ 10"'*, 1 ■ 10~^,3 ■ 10~^ (solid, dashed, dotted curves, 
respectively) Bottom figure is for7 = 7r/6and(5= 1-10"'', 2-10"'*, 3-10"*, (solid, dashed, dotted 
curves, respectively). Chain line in both cases gives the theory II4-8I I with classically computed 
(Tel. Note that the largest 5 = 3- 10~^ case in the top figure (dotted curve) corresponds to 
T ~ 2, so it is already over the upper border of the regime (c, subsect. 12.3^ 5 > (5mix but the 
agreement with the theory 14-81 1 is still quite good, appart from oscillations. This is due to the 
oscillatory nature of time-correlations making the factorization assumption l)2T6|l justified (on 
average) even for much smaller time t as required. 
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Fig. 6. Finite time quantum correlation sum a{t) = '^l7f,i^QC{t' ,t")/2t (solid curves) and the 
corresponding classical sum ad 

(*) = E^rt"=oCci(t',t")/2t (dashed curves saturating at tJci and 
ending at t ~ 1000) for a = 30, J = 1500. Upper curves are for 7 = 7r/6 while lower curves 
are for 7 — n/2. Chain lines are best fitting asymptotic linear functions corresponding to Ct, 
0.0077t/J for 7 = Tv/2 and 0.103t/J for 7 = 7v/6. 

Fpert(t) = exp (-452 JcJdt^). (4-9) 

We numerically computed C (j2-18|) for J = 1500, a = 30 in order to show that it 
is given by the theoretical value ()2-38|) . The quantum correlation function has been 
computed C{t',t") = {Vt'^t") by means of a traceless perturbation V = ^{Jz/J)"^ — 
^[(2J + 1)(J + 1)/J2]1. In fig. El we show a finite time correlation sum a{t) = 
^ ^^r^^^Q C(t', t") which exhibits a crossover, at the Heisenberg time tp = J/2, 
from the plateau given by cJci to a linear increase Ct due to finite size correlation 
average ()2-38|) C = Aa^x/J- The excellent agreement between prediction H4-9|) and 
full numerical calculation of fidelity is shown in fig. [7| In view of our findings this 
so-callecP^ perturbative regime can be understood as a simple consequence of a 
finite Hilbert space dimension. For times larger than the Heisenberg time tp every 
quantum system behaves effectively as an integrable one, i.e. with a finite time 
average correlation plateau. 

Next, we turn to the regime of nonergodic, say regular classical dynamics of the 
kicked top, which is realized for small value of a. If the classical phase space has a 
mixed (KAM) structure, the non-mixing regime of fidelity decay may be obtained 
by choosing a localized initial state (e.g. coherent state) located in a regular part 
of the phase space. Such a situation may easily lead to the opposite conclusion 
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Fig. 7. Quantum fidelity in the perturbative regime S < 5p for a — 30, J — 1500, and 7 = 7r/2 and 
7r/6, calculated as a full trace Hilbert space average as a function of the scaled variable J{St)^. 
For 7 = 7r/2 data for 3=1- 10~® (solid curve) and 5 ■ 10~^ (dotted curve) are shown. For 
7 = tt/6, (5 = 3- 10"^ (solid) and 1 • 10"'^ (dotted) are shown. Note that for both 7 the curves 
for both 5 practically overlap. The chain curves are theoretical predictions 14-91 1 with classically 
computed a^- 



(as compared to generic situation) for an insufficiently large dimension M. As dis- 
cussed in subsect lT^ the fidelity fluctuation plateau is determined by the number 
of constituent propagator eigenstates \4>n) which are effectively needed to expand 
the initial state. For a coherent state sitting inside a (not too large) regular (KAM) 
island this number can be fairly small for numerically realizable Hilbert space di- 
mensions, thus prohibiting any significant fidelity decay as observed in Ref:® We 
would still see the initial quadratic decay in the linear response regime but we would 
not be able to verify higher orders in the long-time expansion of fidelity. In order 
to produce a situation numerically as clean as possible, we choose a small value of 
parameter a = 0.1, such that the classical dynamics is almost integrable and the ma- 
jority of phase space corresponds to regular motion so that the number of constituent 
eigenstates for coherent states is as large as possible (on average). 

Here we focus on the case 7 = 7r/2. For small a, the quantum and classical 
evolution is a (slightly perturbed) rotation around y axis and the time averaged 
perturbation can be computed analytically (to leading order in a) as 

V = -^{J^ + J^) = \{l- (Jy/Jf) , v = \{l-y'). (4-10) 
We will now use these approximate analytical results for a — > to compare with 
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numerics for a = 0.1. Note that our leading order analytical approximations could 
easily be systematically improved using a classical perturbation theory (treating a 
as a perturbing parameter). Since the agreement, as shown below, is almost perfect 
in all cases, we see no need for refinement at this level. 

First consider a random initial state. Starting from expression (|2-21|) . F{t) = 
|(exp {itV6/h))\'^ , write the fidelity as a sum over all eigenvalues of Jy , namely (2m — 
1)^, for m = 1, . . . , J/2 (in OE subspace), 

2 

(4-11) 

For large J we can replace the sum with an integral and get 

2 

(4-12) 

where erfi(z) = JJ^ e~^^ dt is a complex error function with an asymptotic limit 



Fit) 



J/2 



- exp {i6t{2m - l)'^/iJ) 
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Fig. 8. Fidelity in the near-integrable regime for a = 0.1, 7 = n/2, S = 0.01, J = 100, and 
p = 1/J, in the OE subspace. Solid curve gives the result of a numerical simulation. Isolated 
dots denote differences between numerical calculation and the analytic formula 14-1211 for a — > 
|Jiium.(t) — Jknaii.(i)|. The dotted line gives the predicted asymptotic decay oc t^^ , and the 
dashed/chain curves are the predicted fidelity decays at small times, namely the second order 
expansion F{t) = 1 - [JtSf /ISO, and the 'improved' Gaussian result F{t) = exp(-( Jt5)Vl80). 

lim^^-^oo |erfi(ie*'^/^A/x)| = 1 to which it approaches by oscillating around 1. We thus 
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have an analytic expression for the fidehty (|4-12|) in the case of an uniform average 
over the Hilbert space or, equivalently, for a random initial state. Its asymptotic 
decay is t^^ which agrees with the general semiclassical asymptotics ()2-26|) . We 
expect initial quadratic decay ()2-19|) for small times t < T^e- The decay rate r^e 
is determined by the time averaged correlation C ()2-18|) which can be calculated 
explicitly in the limit a — > where the classical correlation function Cci(t) alternates 
for even/odd times as Cci(2t) = {2^{0)z'^{2t))/4 = -1/90, Cci(2t + 1) = {z'^{0)z^{2t + 
l))/4 = 1/45, giving 

Cc,l„„ = 1 (-i; + 1) = _L. (4.13) 

Fidelity is expected to decay as ()2-19|) with Tne = VlSO/ {J6), for short times, t < r^e. 
The short-time formula ()2-19() and the full analytic expression H4-12() are compared 
with the numerical simulation in fig. El The agreement is very good and, surprisingly 
enough, the Gaussian approximation F{t) = exp(— (t/rne)^) for small times is ob- 
served to be valid considerably beyond the second order expansion (|2-19|) . Quite in- 
teresting is the regime where the decay time Tne = \/l80/ (JS) for a "regular" dynam- 
ics with a random initial state will be smaller than the decay time = 1/(5^J^(Tci) 
(|4-8|) for a "chaotic" dynamics. This will happen for 6 < l/( JdciVTSO). This border 
has the same scaling with J as Jmix <|2-43() . 




50 100 150 200 250 300 350 400 450 500 

t 



Fig. 9. Fidelity for 7 = it/2, S — 0.025 and J — 100 on the OF subspace. The dashed curve is a 
simulation for q = 30 (mixing regime, full trace average). The pluses are for a pure coherent 
initial state (see text for details) at a = 0.1 (non- mixing regime). The chain and solid curves 
are, respectively, the theoretical exponential 14-8II and Gaussian 14T6I I decays. 
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Second we consider an SU(2) coherent initial state (|4-5|) . We could perform an 
exact analytical calculation for the fidelity decay in this particular case. Rather 
than performing this calculation, we will illustrate the usefulness of a semiclassical 
formula for (E^TJ). This is a more general approach, as an explicit analytical 

calculation is usually not possible. Let us denote by (p the spherical angular 
coordinates measured with respect to the y-axis. Then (/ = cos-!? = y,<f) represent 
canonical action- angle coordinates for the integrable case a ^ 0. Furthermore, the 
coherent state ()4-5l) acquires a semiclassical Gaussian form (|2-27() in the EBK basis 
|n), Jy\n.) = n\n), namely 

|(n|^,(^)|'ocexp ^ . , fi=- (4-14) 

The squeezing parameter A reads 

A = l/sm^^ = l/{l-y'^). (4-15) 

In order to apply the general formula (|2-31() we need to express the classical time 
average (|4-lUp in terms of a canonical action, v{I) = (1 — /^)/4, and evaluate the 
derivative, \dv{l)/dl\^ = giving = exp {-d^Jt^I^{l - l'^)/8). Rewriting 

this expression in terms of original spherical angles and (/? measured with respect 
to z-axis, we obtain 

= exp i 2^ '^i^^ + ^i^^ ^ '^^^ ^V'i ) • (^'l^) 

We checked this formula by means of a numerical simulation and the result for 
a coherent state centered at (i?,<^) = 7r(l/\/3, l/\/2), for which sin^ 219 sin^ + 
sin^ 1? sin^ 299 = 0.96, is shown in fig. |UJ We can see that for t > 50 the fidelity 
in a non-mixing regime (rne-coh = 23) is lower than the fidelity in a mixing regime 
(Tern = 42). For larger times, t > t*, non-mixing decay F^°^ displays revivals of 
fidelity which are a simple manifestation of a beating phenomenon due to sligtly dif- 
ferent frequencies of revolution of wave-packets along unperturbed and perturbed 
KAM tori. Note that this effect of revivals is drastically reduced in more than one 
dimension, d >2, where there are typically several different incomensurate frequen- 
cies .^^ 

Next we make numerical experiments on the decay of purity fidelity. Again we 
choose the JC Hamiltonian and treat the oscillator and the spin as two parts of a 
composite system. We first report a calculation with a strong perturbation 6 = 0.1, 
that rapidly exceeds the realm of validity of linear response, in Fig. where the 
main figure gives the purity fidelity and the inset the fidelity for comparison. For 
the fidelity decay (inset) we find excellent agreement with exponential decay (|2-17|) 
in a chaotic regime and with a faster Gaussian decay in a regular regime ()2-31|) . 
where the decay rates are fixed as above. For purity fidelity, however, we find 
already at f ~ 20, that the decay starts to be infiuenced by the saturation value 
of Fp(t 00) ~ 1/(2J + 1). Therefore purity fidelity is higher for the integrable 
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Fig. 10. Purity fidelity Fp{t) (main figure) and squared fidelity [F{t)]'^ (inset) in chaotic regime 
(solid curves) and in integrable regime (dotted curves), for 5 — 0.1. The dashed lines indicate 
the linear and quadratic approximation respectively. Note the differences in vertical scales. 



case than for the chaotic one not only at short times, as expected, but even at large 
times. This is relevant, because we shall next choose a weak perturbation 5 = 0.005 
to avoid this problem. We expect and find the crossover after a fairly short time. 
This calculation allows comparison with theory as well as an illustration of the 
evolution of the square of the Wigner function, corresponding to the reduced density 
matrix pi (t) for the angular momentum states on the sphere using the definition of 
Ref.^ Near the top and bottom of Fig. ^2 we see this evolution for the chaotic 
and the integrable Hamiltonian respectively. In the center of the figure we plot the 
purity fidelity on the same time scale as the Wigner functions in the main frame and 
an amplification of short times in the inset. Again, we note faster decay of purity 
fidelity in the case of regular dynamics as compared to chaotic dynamics. We observe 
detailed agreement of numerics with results obtained from the numerical values of 
the correlation integrals ()2-9l3-5|) reproducing the oscillatory structure. From the 
same correlation integrals we obtained the coefficients for the linear and quadratic 
decay, which agree well if we discard the oscillations. We see a crossing of the two 
curves at t = tp ~ 12 for Fp. These times are larger than the Zeno time (~ 1) 
and indicate the competition of the decay rate and the decay shape as expected for 
a non-small value of h = 1/4. It is important to remember that the integral over 
the square of the Wigner function gives the purity and therefore the fading of the 
picture will be indicative of the purity decay. On the other hand the movement of 
the center is an indication of the rapid decay of fidelity (not shown in the figure) . 
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Fig. 11. Echo dynamics for weak coupling: 5 = 0.005. Square of the Wigner function for chaotic 
dynamics (top diagrams) and integrable dynamics as a function of time (bottom diagrams) at 
times corresponding to the axis. Grey code: top left. Purity fidelity is shown in the frame on 
the same time scale and for short times in the inset. Full curves show the complete numerics, 
symbols the evaluation starting from the numerical correlation functions of Fig. |21 and dashed 
curves the linear or quadratic approximation. 



§5. Summary and discussion 

In this paper we have outhned a simple theory of stability of quantum time 
evolution with respect to small variation of the Hamiltonian or, more generally, of 
the time-evolution propagator. The central object of study is the fidelity, a cross 
correllation between two states originating from the same initial state but subject 
to two slightly different time evolutions. An alternative interpretation of the whole 
theoretical setup can be given in terms of quantum dynamical irreversibility in the 
Loschmidt echo gedanken experiment. This intepretation allows even for more gen- 
eral characterization of the quality of quantum echo, as given for example by purity 
fidelity which measures the quality with which the final state, after an echo exper- 
iment, factorizes in the direct product Hilbert space of a composite system. When 
one part of a composite system is interpreted as environment, then the decay of 
purity fidelity can be directly linked to the growth of decoherence. 

The main result reported here is a simple linear response (or Kubo-like) formula 
H2-9|) which relates the fidelity decay, and quite similarly also the purity fidelity 
decay, to the total sum (or integral) of two point time autocorrelation function of 
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the perturbation. In the hmit of infinite Hilbert space dimension we have found 
exponential fidehty decay on a time-scale Tem oc for quantum ergodic and 

mixing systems, whereas for non- ergodic systems we have found much faster decay 
(in the sufficiently 'quantum' regime where 6 <^ h, but note that both, 6 and h are 
dimensionless parameters) on a time-scale r^e (x h/6 for random initial states or 
Tne-coh OC. H^^'^ /5 for Coherent initial states (minimal uncertainty wavepackets) where 
the fidelity F(t) is given by a Fourier transformation of the local density of states 
of the time averaged perturbation operator. A special emphasis has been given to 
the semiclassical theory of fidelity of small but finite values of h, where different 
regimes and the corresponding time and perturbation scales are carefully discussed, 
and where fidelity decay may asymptotically (as ^ — > 0) be evaluated in terms of 
classical quantities only. Interestingly, finite size fiuctuations of fidelity (for very 
long times at a finite Hilbert space dimension) have been shown to be given by the 
inverse participation ratio of the eigenstates of the perturbed evolution operator in 
the eigenbasis of the unperturbed propagator. The surprising aspects of our relations 
are basically due to non-interchangability of the limits 5^0 and h ^ as the 
relevant decay time-scales are only functions of the ratio 6/h. Therefore a different 
and intuitively expected behaviour, namely faster fidelity decay for mixing than 
regular dynamics, is obtained in the 'classical' regime where h <^ d (or making the 
limit ^ — > prior to 6 0). We should note that a similar, reassuring result has been 
found by applying our formalism to inspect the analogous classical fidelit}^'^'^ 
for the unitary Perron- Frobenius evolution of volume (area) preserving mapsP there 
the classical fidelity for regular dynamics has been found to decay on a time-scale 
(X 5~^, i.e. the same as for quantum fidelity decay of regular dynamics, whereas for 
a chaotic dynamics, the fidelity decay is governed by a maximal Lyapunov exponent 
A on a short time-scale oc ln(l/(5)/A. 

Our findings may have some important implication in the understanding of the 
stability of quantum computation^- where different regimes and time-scalesi^ should 
perhaps be explained in terms of intrinsic dynamics of a particular quantum algo- 
rithm. On the other hand, our results may also shed new light on the relation 
between decoherence and dynamics.!^ In particular, due to the quite special role of 
coherent initial states on the short Ehrenfest time-scale one may expect the results 
to be quite different for a random initial state and/or for longer time-scales (se^^). 
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